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Abstract 

If Ultra-High Energy Cosmic Rays (UHECRs) are originated from nearby galaxies, 
modeling of the distribution of nearby galaxies is important to an accurate estimate 
the source number density of UHECRs. We investigate uncertainty of the selection 
function of the Optical Redshift Survey (ORS), which we used to construct a source 
model of UHECRs. The investigation is based on a comparison of numbe counts of 
ORS galaxies with those of the spectroscopic sample of the Sloan Digital Sky Survey 
(SDSS) Early Data Release (EDR). We carefully count galaxies in the same absolute 
magnitude bin from the two samples. We find a slight systematic overestimate of 
the ORS counts outside 5000 km s" 1 by about a factor of 2. We revise the selection 
function of the ORS assuming that the SDSS counts are correct. Our revision is 
based on the absorption given in the ORS catalog as well as that computed from 
Schlegel et al. (1998), which is systematically larger than the former by As ~ 0.1 
mag in the region of low absorption. It is found that introduction of Schlegel et 
al.'s absorption changes one of the parameters of the ORS selection function by more 
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than 10%. The revision should be taken into account in the future analysis of the 
source number density of UHECRs based on the ORS. Using the revised selection 
function, we determine the global structure of the Local Supercluster (LSC) with a 
source model of UHECRs, that is, a number- density model consisting of a uniform 
spherical halo and an exponential disk. We find that the revision is insignificant in 
terms of the structure of the LSC. However, the revised selection function will be 
useful to other studies such as peculiar velocity and correlation function. 

Key words: galaxies: Local Group — galaxies: general — galaxies: structure - 
methods: data analysis — surveys 



1. INTRODUCTION 



Since galaxies are tracers of the underlying matter distribution in the universe, they are 
often used as sources of Ultra High Energy Cosmic Rays (UHECRs), i.e., cosmic rays with 
energy above 10 19 5 eV, in conventional acceleration scenarios (Ide et al. 2001; Blanton et 
al. 2001; Smialkowski et al. 2002; Yoshiguchi et al. 2002). In such scenarios UHECRs are 
considered to be charged particles. Sources of UHECRs must be located in the limited volume 
within a typical radius ~ 50 Mpc because of interaction of such particles with CMB photons 
(GZK effect) (Greisen 1966; Zatsepin and Kuz'min 1966). On the other hand, cosmic rays below 
10 19 ' 5 eV are able to come from much larger volume of radius ~ 100 Mpc, which is determined 
by deflection in the intergalactic magnetic field. The cosmic ray spectrum should therefore 
show a sharp cut off at 10 19 ' 5 eV (GZK cutoff). However, the observed spectrum is featureless 
and the GZK cutoff can not be seen (Takeda et al. 1999). There is another problem concerning 
arrival directions of UHECRs. No significant large scale anisotropy is found in arrival directions 
of 47 UHECRs observed so far. However, one triplet and three doublets with a separation angle 
of 2.5° are observed among the 47 events. Chance probability to observe such clustered events 
under an isotropic distribution is smaller than 1 % (Takeda et al. 1999). 

In order to resolve these problems, there are several groups which perform numerical 
simulation for the energy spectrum and arrival directions of UHECRs. Previous studies can be 
divided roughly into two types. The first uses actual galaxy samples as realistic source models 
of UHECRs as mentioned above (Ide et al. 2001; Blanton et al. 2001; Smialkowski et al. 2002; 
Yoshiguchi et al. 2002). The second adopts simple source models which represent the Local 
Supercluster (LSC), such as a two dimensional sheet (Sigl et al. 1999) or a highly flattened 
oblate spheroid (Lemoine et al. 1999). 

A major difficulty in calculating the UHECRs flux numerically with an actual galaxy 
sample is that galaxies in a given magnitude-limited sample are biased tracers of matter dis- 
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tribution because of the flux-limit. In our previous study, we examined whether current obser- 
vations of UHECRs could be explained by conventional acceleration scenarios (Ide et al. 2001; 
Yoshiguchi et al. 2002) using a whole sky survey, Optical Redshift Survey (ORS; Santiago et 
al. 1995). However, the correction for the selection effect is often very large. For example, the 
selection function of the ORS given in Santiago et al.(1996) is roughly equal to 0.01 at 8000 km 
s . Blanton et al.(2001) showed that the small local overdensity which they derived from the 
Center for Astrophysics Redshift Survey and the IRAS PSCz Survey gave the predicted GZK 
cutoff which is discrepant with the current observation at only 2 o level. As is evident from 
this result, it is important to know uncertainty of the large correction. 

A major concern with studies of the second type is that their source models are too 
primitive. It should be emphasized that UHECRs above 10 19,5 eV must originate within about 
50 Mpc from us. Accordingly, the spatial distribution of nearby galaxies plays an important 
role for the distribution of arrival directions of UHECRs. However, in addition to the primitive 
form of their models, parameters such as radial and vertical length of the models are not well 
constrained. 

In this paper, we investigate the uncertainty of the ORS selection function by comparing 
the galaxy number counts of galaxies as a function of recession velocity between Sloan Digital 
Sky Survey: Early Data Release (SDSS EDR) (Stoughton et al. 2002) and the ORS. The 
SDSS is the imaging and spectroscopic survey in five optical frequency bands with a limiting 
magnitude which is much deeper than that of the ORS, but covering a very limited area of the 
ORS. 

The plan of this paper is as follows: We introduce the SDSS and ORS data in Section 
2. Effect of Galactic absorption is described in Section 3. Comparison of the galaxy number 
counts and revised parameter of the ORS selection function are given in Section 4. In Section 
5, we describe a source model of UHECRs, that is, a model of the galaxy number density, the 
method of determining its parameters, and resulting parameters. In Section 6, we summarize 
the results and discuss their implications. 

2. THE GALAXY SAMPLES 

2.1. THE ORS SAMPLE 

The ORS is drawn from three catalogs, the European Southern Observatory Galaxy 
Catalogue (ESO) covering the southern sky (5 < —17.5°), the Uppsala Galaxy Catalogue (UGC) 
in the northern sky (5 > —2.5°), and the Extension to the Southern Galaxy Catalogue (ESGC) 
in the strip between the ESO and UGC region (—17.5° < 5 < —2.5°). The last catalogue 
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contains no magnitudes other than very poor ones from eye inspection. The sample selection 
and galaxy distribution are described in detail by Santiago et al. (1995). This survey contains 
two subcatalogs, one complete to a B magnitude of 14.5 (ORS-m) and the other complete to a 
B major axis diameter of 1/9 (ORS-d). In this study we calibrate the ORS selection function 
using the SDSS EDR which is a magnitude-limited survey. Accordingly, we use the magnitude- 
limited subcatalog, i.e., ORS-m, which contain 5716 galaxies. The subcatalog ORS-m itself 
is also composed of two subsamples, ESO-m (2437 galaxies) and UGC-m (3279 galaxies). We 
further restrict ourselves to galaxies within 8000 km s _1 and with Mb — 5\og w h < — 0.5 — 51ogiof , 
which corresponds to the apparent magnitude limit of m# = 14.5, where h is the dimensionless 
Hubble constant h = Hq/(100 km s -1 Mpc -1 ) and v is the recession velocity. Very faint nearby 
galaxies with Mb > —14.0 are excluded. Galaxies in the zone of avoidance (|6| < 20°) are also 
excluded. Our final ORS sample contains 5178 galaxies with 2346 galaxies in the ESO-m and 
2832 galaxies in the UGC-m. 

Figure 1 shows the sky distribution of our ORS galaxies in the equatorial coordi- 
nate. Figure 2 shows the absolute magnitude plotted as a function of the recession velocity. 
Throughout the paper, we assume that there are no departures from uniform Hubble expansion 
for simplicity. We use h = 1 when we do not specify h explicitly. An alternative value for h 
changes all absolute magnitudes by 51ogio/i. 

In order to compare the galaxy number counts of the ORS with those of the SDSS EDR, 
a selection function has to be introduced. This selection function of the ORS is derived for 
several subsamples in Santiago et al.(1996) using galaxies to 8000 km s _1 . They assumed that 
the selection function is expressed as 

0ORs(r) = ( fe) feSf) r>rs , (1) 

[1 r <r s 

where r represents distance from our Galaxy, and they set r s to be 500 km s -1 . Values of the 
parameters determined by Santiago et al.(1996) are a = 0.40, b = 6.25, r* = 11100 (km s _1 ) 
for ESO-m and a = 0.36, b = 9.38, r* = 14060 (km s" 1 ) for UGC-m. The flux limit m B = 14.5 
corresponds to Mb = —14 at v — 500 km s _1 , where the selection function is unity. In Figure 3 
we show the selection function of the ORS given by Santiago et al.(1996) as a function of 
recession velocity together with that of the SDSS (explained in the next subsection). 

2.2. THE SDSS SAMPLE 
2.2.1. Completeness 



The SDSS is the imaging survey in five bands (u',g',r',i',z') followed by spectroscopy. 
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The spectroscopic sample of the SDSS EDR contains 1504 galaxies with measured recession 
velocity v < 8000 km s" 1 . The region covered by the spectroscopy of the SDSS EDR is shown 
in Figure 4. In the following, we refer to the region at § > 50° as region 1, two regions in the 
equatorial stripe as regions 2 and 3. All the data in these three regions are collectively referred 
to as the ALL region. 

At the brightest end, it is suspected that problems in the deblender and/or possible 
saturation limit the completeness of the SDSS. There are 130 ORS galaxies which fall in the 
ALL region. We tried to match these 130 galaxies with the SDSS catalogs. We found that 72 
ORS galaxies out of the 130 could be matched with the SDSS spectroscopic objects and that 
43 could be matched with the photometric objects (but not with the spectroscopic objects). 
Remaining 15 galaxies are not included either in the spectroscopic catalog nor in the photo- 
metric catalog. Among them, 10 suffer from problems in the deblender, 1 is too bright (~ 9 
mag), and remaining 4 galaxies are missed by some unknown reason. Thus the completeness 
of the SDSS spectroscopic catalog is quite low, 72/130 = 55%, at the ORS limiting magni- 
tude. However, if supplemented by the photometric catalog, it increases to 115/130 = 88%. 
The SDSS spectroscopic catalog have a very high completeness for galaxies including very faint 
ones (above 99% for r' < 18) (Strauss et al. 2002). We add the 43 SDSS photometric objects 
in the above 1504 galaxies assigning them the recession velocities given in the ORS. Thus, our 
SDSS sample now includes 1547 galaxies. We extract our final sample from these 1547 galaxies 
to be compared with the ORS sample including the 130 galaxies. 

2.2.2. Conversion from tub tog' magnitude 

The ORS B magnitude is based on the photographic plates and its bandpass is similar to 
that of Couch-Newell Bj (Couch and Newell 1980). We use only SDSS g'-band magnitude in this 
study because it is close to the Bj magnitude. In order to carefully count galaxies in the same 
magnitude range, we have to know magnitude difference between m# and g' . For this purpose, 
we made a direct comparison of the magnitudes of the matched 72 plus 43 galaxies between 
the two surveys. However, variance of the magnitude offset was very large (g' — nis ~ 1.5 at 
the maximum). At such bright end (the ORS-m is complete to ms = 14.5 mag), it is suspected 
that problems in the deblender limit the accuracy of the SDSS magnitude. Accordingly, we 
derive the magnitude offset in the method explained below. 

There are three differences of the magnitude system between the SDSS and the ORS. 
First, response functions of these two photometric band systems differ from each other. 
Magnitude offset due to this difference can be computed from galaxy colors in various pho- 
tometric band systems given in Fukugita et al.(1995) for typical morphology types grouped 
into E, SO, Sa-Sb, Sb-Sc, Sc-Sd, and Sdm-Im. The ORS gives the morphological type index 
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to individual galaxy, whereas the SDSS does not. However, Strateva et al.(2001) showed that 
galaxies have a bimodal u* — r* color distribution corresponding to early (E, SO, Sa) and late 
(Sb, Sc, Irr) morphological types and that the two types can be clearly separated by a u* — r* 
color cut of 2.22 independent of magnitude. Accordingly, we compute the magnitude offset 
for early and late morphological types and use them throughout the paper. Averaging over 
morphology types of galaxies in our ORS sample, we adopt g' — Bj = —0.35 and —0.24 for early 
and late type galaxies, respectively. 

Second, there is a difference about the method for measuring the flux. In order to 
measure a constant fraction of the total light, independent of Galactic absorption and cosmo- 
logical dimming, the SDSS has adopted a modified form of the Petrosian (1976) magnitude. 
Magnitudes are measured within the aperture of twice the Petrosian radius, which is insensitive 
to absorption or cosmological dimming. The detailed explanations are presented in Blanton 
et al.(2001), Yasuda et al.(2001), and Stoughton et al.(2002). On the other hand, the ORS is 
based on an isophotal magnitude. The relationship between the SDSS Petrosian magnitudes 
and isophotal magnitudes is examined in Blanton et al.(2001). Difference between the two 
magnitudes Am = m petro — mi SO ranges from —0.2 to 0.1 according to several choices of the 
isophotal limit in case of r' band. It is, however, not known exactly what value we should take 
for Am in the g' band. Here we tentatively take Am = —0.1 with an estimated uncertainty of 
±0.2. 

Finally, SDSS magnitudes are not the traditional logarithmic magnitude but inverse 
hyperbolic sine magnitude, which is described in detail by Lupton et al.(1999). However, these 
magnitudes differ by less than 1 % in flux above g' = 22.60. We neglect this difference throughout 
the paper. 

Taking account of these three differences between the SDSS and the ORS magnitudes, 
we take the total amount of the magnitude offset as g' — Bj + Am = —0.45 ±0.2 and —0.34 ±0.2 
for early and late type galaxies. 

2.2.3. Selection of galaxies from our SDSS sample 

We apply the following selection criteria to our 1547 SDSS galaxies with v < 8000 km 
s -1 . The limiting magnitude is set at g' = 17.65, which was adopted in the luminosity function 
study by Blanton et al.(2001). This corresponds to the limit in the absolute magnitude of 
M g i = 2.65 — 51ogiof • The limiting magnitude of the ORS is Mb = —14.0, which corresponds to 
M g i = —14.45 and —14.34 for early and late type SDSS galaxies. Accordingly we also exclude 
galaxies fainter than this limiting magnitude. The absolute g' band magnitude plotted as a 
function of the recession velocity is shown in Figure 5. The sudden change in the density of 
galaxies which occurs at around 5000 km s _1 reflects a real large scale structure. Among the 
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1547 galaxies we find a tiny fraction (~ 3%) of galaxies with g' 2> 17.65, which may be assigned 
unrealistically faint absolute magnitudes because of some problem in deblending. There is 
concern that these galaxies are actually more luminous than g' = 17.65. They are not included 
either in the sample or in Figure 5 because of the above cut in the absolute magnitude. However, 
rejection or inclusion of this tiny fraction of erroneous data does not affect our results. Our 
final SDSS sample contains 1078 galaxies. Figure 6 shows the distribution of galaxies in our 
SDSS sample and ORS galaxies in the EDR region in the equatorial coordinate. The numbers 
of observed galaxies in our SDSS and ORS sample within each region are tabulated in Table 1. 
Note that the whole region covered by the SDSS sample is located in the area covered by 
UGC-m (See Figure 1). 



2.2.4- Selection function 

The selection function of the SDSS can be computed by using the luminosity function 
(Blanton et al. 2001). The luminosity function can be expressed by the Schechter function as 

$(M) oc Kr°- 4 (i*-M.)(«+i) x exp [i(r - 4 ( A/ - M *)] , (2) 

where M is the absolute magnitude. The parameters for g' band determined by Blanton et 
al.(2001) are M* = -20.04 ±0.04 mag, a = -1.26 ±0.05 for tt m = 0.3 , fi A = 0.7 cosmology. 
When we restrict our analysis to galaxies with absolute magnitudes M < Mo, the selection 
function is given by 

J_^$(M)dM 

where M max (r) =min(M ,Mi im (r)) and Mi im (r) is the limiting magnitude for a galaxy at distance 
r to be included into the sample. We set Mq to be —14.45 and —14.34 for early and late type 
galaxies which are the absolute magnitudes in g' band corresponding to Mb = —14. The SDSS 
selection function obtained in this way is shown in Figure 3. 



0sdss(h= m,^,^,,, , (3) 



3. GALACTIC ABSORPTION 



Galactic absorption introduces a nonuniform, direction-dependent selection effect. The 
effective magnitude cutoff limit for a galaxy lying in the direction (/, b) is given by m\\ m = 
^iim,obs — A\(l,b), where A x is the absorption at a given band and (l,b) is the Galactic coordinate. 
The selection function with the absorption taken into account, is now given for each of the 
surveys as 
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4>obs(r,l,b) 



L 



max(r,I,6) $ ( M ) dM 



(4) 



J^$(M)dM 



where 



M max (r,l,b) 



M max (r)+A x (l,b) 
M max (rx 10 a2A ^). 



(5) 
(6) 



Note that M max (r) is related to the flux limit /i; m as M max (r) = const — 2.51og 10 47rr 2 /i im . Thus 
we can obtain the selection function o b s (r, /,6) in which the absorption is taken into account 
from the one 0(r) which incorporates no absorption as 



The ORS gives the absorption Ab to individual galaxy based on Burstein and Heiles 
(1982) while the SDSS gives the absorption A g > based on Schlegel et al.(1998). As shown in 
Fig 7, a comparison between the ORS Ab and the SDSS A g i reveals a systematic difference 
in the region of low absorption where Ab ~ 0.0 while A g > ~ 0.1. Accordingly we make two 
computations, one with Ab given in ORS and the other with Ab computed from Schlegel et 
al.(1998) as Ab = 4.1E(B — V), in the comparison of the galaxy number counts. The absorption 
for the SDSS galaxies are computed by A g > = 3.793E(B — V) (Stoughton et al. 2002). 

The ORS is based on the isophotal magnitude as mentioned above. Absorption dims the 
isophotal magnitude because not only every point in the galaxy become dimmer, but the radius 
of the isophote shrinks. In Santiago et al.(1996), considering these effects, observed magnitude 
is given by m B obs = m B +'y m A B (l,b), where 7 m is a free parameter of the selection function to 
be determined from the data themselves. We expect 7 m > 1.0 in the isophotal magnitude, the 
exact value depending on the two-dimensional light profile of the object. (If m B obs is the total 
magnitude, 7 m = 1.0.) The parameters 7 m derived in Santiago et al.(1996) are 1.0 for ESO-m, 
0.56 for UGC-m. However Santiago et al.(1996) noted that errors for 7 m are large and these 
values are consistent with j m — l. Since we can not adopt a value of 7 m greater than unity 
appropriately, we set j m — 1 for the ORS throughout the paper. The SDSS adopt Petrosian 
magnitude as mentioned in Section 2.2. In this case, the Petrosian radius of the object does 
not shrink. Thus we also set j m = 1 for the SDSS throughout the paper. 

4. CALIBRATION OF THE ORS SELECTION FUNCTION 

In this section, we compare the galaxy number counts derived from our SDSS sample 
(1078 galaxies) with those from our ORS sample (130 galaxies). We start by counting the 
number of the observed galaxies in the SDSS and in the ORS, which lie in each of the regions 
covered by the SDSS, in velocity bins of width 1000 km s _1 . Then, we correct this number for the 



(j) ohs (r,l,b) = 0(r x 10 



,0.2A x (l,b) 



(7) 
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selection effect using respective selection functions shown in Figure 3. The galactic absorption 
for the ORS is computed following Schlegel et al.(1998). The results are summarized in Table 1 
and shown in the upper panels of Figure 8. The error bars are defined so as to represent a 
square root of the number of galaxies in each velocity bin. Since observed galaxies in the SDSS 
are much more than in the ORS, we show errorbars only for the ORS. We find that number 
counts of the ORS are systematically larger outside 5000 km s" 1 . 

Cosmic variance is always a concern in this kind of analysis. For this reason, we counted 
the galaxies separately in each of the three regions. As shown in Table 1 and Figure 8, we 
found the overestimate of the ORS in all the three regions. Although we can not draw strong 
conclusion as for the significance of this overestimate considering the current small area coverage 
of the SDSS, accumulation of the data may allow us to do more extensive investigation in the 
near future. For the moment, we regard the systematic difference as significant because we have 
carefully counted galaxies in the same magnitude range (Mb < — 14) and the SDSS spectroscopic 
catalog has a very high completeness except for the bright end as shown in Strauss et al.(2002). 

We revise the selection function of UGC-m by changing the parameter b, which is sensi- 
tive at large distances, so that the ORS galaxy counts coincide with the SDSS counts. The effect 
of changing b is clearly visible in Figure 9, where the galaxy number counts of the ORS in ALL 
region are shown for various choices of the parameter b. The revised values of b for regions 1,2, 
3, and ALL region for the adopted values of the magnitude offset of g' — Bj = 0.34 mag (early 
type) and —0.45 mag (late type) are summarized in the top line of Table 2. Absorption aver- 
aged over ALL region is also given in this table. Comparison of the galaxy number counts using 
the revised selection function are shown in the lower panels of Figure 8. We show the revised 
selection function for ALL region in Figure 10, together with that by Santiago et al.(1996). 

We examine the effect of uncertainty of the offset. The values of b for Am = —0.1 ±0.2 
mag are also shown in Table 2. The b values changes by ~ ±0.3 for each region. Furthermore, 
we also tabulate results derived by using the original absorption given in the ORS catalog 
in Table 2. Replacing the original ORS absorption with that based on Schlegel et al.(1998) 
changes the b values by ~ 0.9, which is larger than the change caused by the different offset 
values. In all cases, the values of b become smaller than that derived by Santiago et al.(1996) 
(W = 9.38). 

5. SPATIAL STRUCTURE OF THE LOCAL SUPERCLUSTER 

It is long known that there is an excess of galaxies in the vicinity of the north galactic 
pole, which is now referred to as the LSC (e.g., de Vaucouleurs 1956; Yahil et al. 1980; Tully 
1982). Tully (1982) illustrated the three-dimensional distribution of galaxies in the LSC and 
described its spatial structure as consisting of two components, one is a flattened disk and the 
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other is a roughly spherical halo. 

Only luminous galaxies are often counted as sources of UHECRs, and faint galaxies are 
neglected. However, the boundary between luminous and faint galaxies is not established at 
all in the context of acceleration scenarios. If luminous galaxies in the LSC are distributed 
differently from fainter ones, the distribution of arrival directions of UHECRs may change 
with the limiting absolute magnitude of the sample considered. Thus we determine the spatial 
structure of the LSC using galaxies of different absolute magnitudes. In Figure 11, we show 
distribution of galaxies more luminous than Mb = —14,-16,-18 in the LSC in the super- 
galactic coordinate. The two components, the flattened disk and the roughly spherical halo, 
are clearly visible. Here we restrict our analysis to galaxies within 2000 km s -1 from the Virgo 
cluster (SGX=SGZ=0 km s -1 , SGY=1000 km s" 1 ). 

We assume that the number density of the halo is uniform while that of the disk is 
exponential both in the radial and the z-directions. Thus, the number density of galaxies of 
our model is expressed as 



where p, z are the cylindrical coordinates in which the Virgo Cluster is located at the origin 
(p = 0, z = 0) and the plane of the disk component coincides with the z = plane. The parameter 
rii in Eq.( 8) represents the ratio of the halo component to the disk component. The parameters 
of a and (3 are the radial and vertical scale lengths of the number density of galaxies in the 
disk component, respectively. 

The probability fi/rt) that a galaxy is observed at position f\ is proportional to the 
number density multiplied by the selection function at ~rt and divided by the integral of these 
values over the region which we consider, that is, the region within 2000 km s _1 (20 h~ l Mpc) 
from the Virgo cluster. Thus we can write 



This probability function is independent of the normalization of the number density. For 
computational convenience we define a likelihood function to be the logarithm of the product 
of Eq.( 9) over all galaxies in the sample and in the region considered: 



We determine the most probable values of a, ft and n\ by maximizing Eq.( 10) with respect to 
these parameters for each of the three samples with Mb < —14,-16 and —18. In the calculation 
of Eq.( 10), we take account of the absorption from Schlegel et al.(1998). The ORS excluded 
the region of \b\ < 20°, and in addition, the region at higher latitude in which the B-band 
extinction according to Burstein and Heiles (1982) is greater than Ab = 0.7 mag (this excludes 
only 0.18 sr). However, a tiny fraction (~ 4%) of our ORS sample lie at the position where the 




(8) 




(9) 




(10) 
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absorption from Schlegel et al.(1998) is larger than Ab = 0.7 mag. We exclude those galaxies 
in the calculation of Eq.( 10). 

In order to determine the most probable values, we need to obtain the revised selection 
function for ESO-m as well as for UGC-m. The revised b value for ESO-m is derived as 
follows. We compute the ratio of the value of the original UGC-m selection function to that 
of the revised selection function at 8,000 km s _1 , and then, calculate the revised b value for 
the ESO-m selection function so that the above ratio takes the same value as for the UGC-m 
selection function. Because we want to determine the most probable values for galaxies with 
different absolute luminosities, we compute the selection function corresponding to each limiting 
magnitude as follows. First, we obtain the luminosity function of the ORS by differentiating 
Eq.(l) with respect to r. Note that the derivative of the selection function gives the luminosity 
function as seen, for example, Eq.(3). With this luminosity function, we calculate the ORS 
selection function by the equation similar to Eq.(3) with various values of Mq. We choose 
Mq = — 16 and —18 as well as the original choice of M = — 14. In Figure 12, the original and 
the revised selection functions of UGC-m and ESO-m samples are shown for Mb < —14,-16 
and —18. It is noted that for Mb < —18 both UGC-m and ESO-m samples are complete to ~ 
3000 km s^ 1 (see also Figure 2). 

We show the values of the likelihood function with respect to a,/3 and n\ in Figures 13. 
The most probable values of the model parameters are summarized in Table 3. It is found that 
the revision of the selection function does not affect the results in this analysis. Most probable 
values of (3 are almost the same for all the luminosity ranges, while a and n\ are larger for 
Mb < — 18 than for Mb < — 14 and —16. In other words, if we do not see dwarf galaxies 
[Mb > —16), the disk looks more extended outward and halo looks denser. These properties 
are also seen in Figure 11. 

Let us calculate a quantitative value of the local overdensity due to the disk component 
in the LSC. We define the local overdensity 5 as the ratio of number density of galaxies averaged 
over the disk component to that of the halo component, which is assumed to be uniform in this 
paper. (We tacitly assume all the galaxies to have same mass.) We define the volume of the 
disk component over which the number density is averaged by 

9 1 <-ln{nx). (11) 



a^Mpc) ^(/i^Mpc) 
The volume corresponds to the region where the first term in Eq.( 8) which expresses the disk 
is larger than the second term which expresses the halo. The local overdensities for different 
luminosity ranges are also given in Table 3. Baker et al.(1998) calculated the density field 
contours on several radial shells with the Gaussian smoothing, and obtained local overdensity 
~ 10 for the Virgo cluster using the ORS galaxies. Strictly speaking, local overdensity derived 
by Baker et al.(1998) is not defined in the same manner as ours. However, the amount of 
the local overdensity derived by our number- density model is consistent with that obtained by 
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Baker et al.(1998). 



6. SUMMARY AND DISCUSSION 



In this paper, we investigate the uncertainty of the selection function of the ORS, which 
we used as sources of UHECRs, by comparing the galaxy number counts of ORS galaxies with 
those of SDSS EDR galaxies as a function of recession velocity We carefully count galaxies 
from the two data samples so that these galaxies are in the same absolute luminosity range. In 
all the regions covered by the SDSS EDR, we find that the ORS predicts systematically more 
galaxies than the SDSS by a factor of ~ 2 outside 5000 km s" 1 . We revise the selection function 
of the ORS so that the predictions of the two surveys become consistent. It is also found that 
replacing the original absorption given in the ORS catalog with that computed from Schlegel 
et al.(1998) changes the b parameter of the ORS selection function by more than 10 %. The 
revision should be taken into account in the future analysis of the source number density of 
UHECRs based on the ORS. 

Using this revised selection function of the ORS, we determine the spatial structure of 
the LSC with a model consisting of a uniform spherical halo and an exponential disk. The 
model parameters are derived for galaxies in three ranges of different absolute luminosities, 
Mb < —14, Mb < —16 and Mb < —18. The most probable values of the scale length of the 
disk in z direction, /3, are almost the same for the three luminosity ranges while the radial 
scale length of the disk, a, and the density of the halo, ni, are larger for Mb < —18 than 
for Mb < —14 and —16. This implies that giant galaxies (Mb < —18) in the disk are more 
extended radially than dwarf galaxies (Mb > —16) and that giant galaxies show a smaller local 
overdensity than dwarf galaxies. 

Let us consider the effect of this difference of the distribution of galaxies in the LSC on 
the distribution of UHECR arrival directions. Observationally, it is reported that there is no 
significant large-scale anisotropy in the arrival direction distribution of UHECRs (Takeda et al. 
1999). If all the galaxies, both giants and dwarfs, are sources of UHECRs, arrival directions 
could be more anisotropic than the case that only luminous giant galaxies are responsible for 
UHECRs. The observed lack of anisotropy may favor the latter case. Provided that UHECRs 
are generated in star forming regions, their sources are expected to trace the distribution of 
late-type galaxies rather than early-type galaxies. It is known that there is strong spatial 
segregation of the Hubble types in the Virgo cluster (Binggeli et al. 1987). Early-type galaxies 
are more concentrated towards the cluster center than late types. Distribution of galaxies of 
different morphological types in the LSC is also an important factor to be taken into account 
in the analysis of arrival direction of UHECRs. 

We restricted ourselves in this paper to the study of galaxies as the sources of UHECRs. 
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The revision of the ORS selection function made in this study is found to be insignificant in the 
determination of the global structure of the LSC The revised selection function obtained in this 
study will, however, be useful to other studies such as peculiar velocity and correlation function. 
The SDSS will produce imaging and spectroscopic surveys in five bands (u',g',r',i',z') over tt 
steragians in the northern Galactic cap, and will collect spectra of approximately 1,000,000 
galaxies (down to r' lim ~ 17.65). We have used in this paper the SDSS EDR which covers only 
less than 1% of all sky. As data obtained by the SDSS accumulate, we will be able to revise 
the selection function of several surveys with better accuracy, and also to determine the spatial 
structure of the LSC using only the SDSS galaxies. 

Acknowledgments 

The authors are grateful to Dr. K. Shimasaku and SDSS group in University of Tokyo for 
useful discussions. This research was supported in part by Giants-in-Aid for Scientific Research 
provided by the Ministry of Education, Science and Culture of Japan through Research Grant 
No.S14102004. 

References 

Baker J.E., Davis M., Strauss M.A., Lahav O., Santiago B.X. 1998, ApJ, 508, 6 

Binggeli B., Tammann G.A., Sandage A. 1987, AJ, 94, 251 

Blanton M., Blasi P., Olinto A.V. 2001, Astropart. Phys., 15, 275 

Blanton M.R., et al. 2001, AJ, 121, 2358 

Burstein D., and Heiles C. 1982, AJ, 87, 1165 

Couch, W.J., and Newell, E.B. 1980, PASP, 92, 746 

de Vaucouleurs, G. 1956, Vistas in Astronomy, 2, 1584 

Fukugita M., Shimasaku K., Ichikawa T. 1995, PASP, 107, 945 

Greisen K. 1966, Phys. Rev. Lett., 16, 748 

Ide Y., Nagataki S., Tsubaki S., Yoshiguchi H., Sato K. 2001, PASJ, 53, 1153 
Lemoine, M., Sigl, G., and Biermann, P. 1999, astro-ph/9903124 
Lupton R. H., Gunn J. E., and Szalay A. S. 1999, AJ, 118, 1406 
Petrosian V., 1976, ApJ, 209, LI 

Santiago B.X., Strauss M.A., Lahav O., Davis M., Dressier A., Huchra J. P. 1995, ApJ, 446, 
457 

Santiago B.X., Strauss M.A., Lahav O., Davis M., Dressier A., Huchra J. P. 1996, ApJ, 461, 38 

Schlegel D., Finkbeiner D., and Davis M. 1998, ApJ, 500, 525 

Sigl, G., Lemoine, M., and Biermann, P. 1999, Astropart. Phys., 10, 141 

Smialkowski, A., Giller, M., and Michalak, W. 2002, astro-ph/0203337 

Stoughton C, et al. 2002, A J, 123, 485 



13 



Strateva L, et al. 2001, A J, 122, 1861 
Strauss M. A., et al. 2002, A J, submitted 
Takeda M., et al. 1999, ApJ, 522, 225 
Tully R.B. 1982, ApJ, 257, 389 

Yahil, A., Sandage A., and Tammann G.A. 1980, Astrophys. J., 242, 448 
Yasuda N., Fukugita M., Narayanan V. et al. 2001, AJ, 122, 1104 
Yoshiguchi H., Nagataki S., Tsubaki S., Sato K. 2002, astro-ph/0210132 
Zatsepin G.T., Kuz'min V.A. 1966, J. Exp. Theor. Phys. Lett., 4, 78 



14 





Region 1 


Region 2 


Region 3 


V 


ORS 


SDSS 


ORS 


SDSS 


ORS 


SDSS 


(km s -1 ) 


N bs N COI 


N ohs 


iVcor 


N ohs 


iVcor 


N oba 


iV cor 


N ohs 


iV cor 


N ohs 


iV cor 


0-1000 


n 







0.0 


2 


3.3 


1 


1.0 





0.0 


1 


1.0 


1000-2000 


n 








0.0 


23 


62.6 


46 


46.0 


8 


20.5 


14 


14.0 


2000-3000 


2 


11.5 


2 


2.3 


6 


26.3 


16 


16.5 


4 


20.9 


26 


28.4 


3000-4000 


8 


59.7 


30 


39.4 


5 


46.1 


23 


31.9 


1 


7.9 


10 


14.0 


4000-5000 


2 


29.9 


10 


17.6 


3 


45.3 


18 


30.5 


6 


81.9 


41 


68.4 


5000-6000 


1 


18.1 


15 


30.3 


14 


398.1 


79 


161.6 


16 


405.5 


177 


353.4 


6000-7000 


2 


91.6 


19 


44.4 


7 


360.7 


124 


297.5 


15 


944.0 


153 


381.5 


7000-8000 


1 


143.2 


36 


103.1 


1 


101.7 


127 


363.3 


3 


335.2 


110 


318.6 


total 


16 


354.0 


112 


237.1 


61 


1044.2 


434 


948.4 


53 


1815.9 


532 


1179.2 



Table 1. Number of galaxies in our ORS and SDSS samples in each region defined in Figure 4. A^bs 
and N cor represent the number of observed galaxies and that corrected for the selection effect using the 
selection functions, respectively. The magnitude offset is taken as g 1 — Bj + Am = —0.45 and —0.34 for 
early and late type SDSS galaxies. 
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Table 2. Parameters b of the ORS selection function for UGC-m revised in this study according to several 
choices of the magnitude offset between the SDSS and ORS and of the absorption effect. The original 
value of b derived by Santiago et al.(1996) is 9.38. The absorptions averaged over all the galaxies in the 
ALL region are also tabulated. Absorption of the SDSS is computed following Schlegel ct al.(1998). 
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Table 3. Parameters of our number-density model of the Local Supcrclustcr calculated using the ORS 
sample of different luminosity ranges within 2000 km s _1 from the Virgo cluster. The local overdensity 
due to the disk component is also shown (see the text). 



Fig. 1. Distribution of galaxies in the ESO-m and UGC-m samples used in the present study (v < 8000 
km s" 1 , M B < -14.0, m B < 14.5, \b\ > 20°). 
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Fig. 2. Absolute magnitude plotted as a function of the recession velocity for galaxies in our ORS sample. 
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Fig. 3. Selection functions of the ESO-m and the UGC-m given in Santiago et al.(1995) and of the SDSS 
(g' band). The SDSS selection function is calculated from the luminosity function given in Blanton et 
al.(2001). 




17 




Fig. 5. Absolute g' band magnitude plotted as a function of the recession velocity for galaxies in our 
SDSS sample. 
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Fig. 6. Distribution of 1078 galaxies in our SDSS sample (left panels) and 130 galaxies in our ORS 
sample (right panels). 
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Fig. 7. Absorptions given to individual ORS galaxies in the regions of SDSS EDR (left panel) and those 
given to SDSS galaxies (right panel) as a function of right ascension. A systematic difference between the 
two in the region of low absorption is clearly visible. 
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Fig. 8. Number of galaxies in velocity bins of width 1000 km s _1 in each region. Numbers from the SDSS 
are corrected using the selection function derived in this study. Numbers from the ORS are corrected by 
the selection function given in Santiago et al.(1995) (upper panels) and by the selection function revised 
in the present study (lower panels). We show error bars on numbers of the ORS galaxies only, which are 
much larger than those of the deeply sampled SDSS galaxies. 
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Fig. 9. Number of ORS galaxies, which fall into the ALL region and are corrected by the selection 
function given in Santiago et al.(1995), for various choices of the parameter b. It is clearly visible that the 
number counts at large distance are sensitive to this parameter. 




Fig. 10. The selection function of the UGC-m with parameters given in Santiago et al.(1995) (dotted 
line) and revised in this analysis (solid line). The revised b value from the analysis in ALL region is used. 
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Fig. 11. Galaxies more luminous than Mb = —14,-16,-18 in ESO-m and UGC-m data sample within 
2000 km s" 1 from the Virgo cluster, which is the dense region near SGZ=0 km s 1 , SGY=1000 km s 1 . 
Our galaxy is located at the origin of the coordinate system. 
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Fig. 12. The selection function of UGC-m and ESO-m samples for M B < -14, -16,-18. As for ESO-m 
data sample, the revised b parameter is calculated so that the ratio of the value of the ESO-m selection 
function with the revised b value at 8000 km s" 1 to that with the original b value given in Santiago et 
al.(1995) becomes equal to the ratio of the UGC-m selection function. 
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Fig. 13. Values of the likelihood function with respect to a (left panels), j3 (center) and n\ (right) for 
luminosity range Mb < — 14 (upper panels), Mb < — 16 (middle), Mb < — 18 (lower). The solid curves 
represents the results based on the revised selection function while the dotted curves are based on the 
original selection function. Absolute values of the likelihood function are normalized arbitrarily. Note that 
for Mb = — 18 both ESO-m and UGC-m are complete samples in the limited- volume of radius ~ 2000 km 
s _1 from the Virgo cluster (See Figure 2). 
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